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D^' Abstract 

W). The Davies-Fulhng model describes the scattering of a massless field by a moving mirror 

in 1 + 1 dimensions. When the mirror travels under uniform acceleration, one encounters 
^ I severe problems which are due to the infinite blue shift effects associated with the horizons. 

I On one hand, the Bogoliubov coefficients are ill-defined and the total energy emitted 

diverges. On the other hand, the instantaneous mean fiux vanishes. To obtained well- 
defined expressions we introduce an alternative model based on an action principle. The 
usefulness of this model is to allow to switch on and off the interaction at asymptotically 
large times. By an appropriate choice of the switching function, we obtain analytical 
expressions for the scattering amplitudes and the fluxes emitted by the mirror. When 
the coupling is constant, we recover the vanishing flux. However it is now followed by 
transients which inevitably become singular when the switching off is performed at late 
time. Our analysis reveals that the scattering amplitudes (and the Bogoliubov coefficients) 
should be seen as distributions and not as mere functions. Moreover, our regularized 
amplitudes can be put in a one to one correspondence with the transition amphtudes of 
an accelerated detector, thereby unifying the physics of uniformly accelerated systems. 
In a forthcoming article, we shall use our scattering amplitudes to analyze the quantum 
correlations amongst emitted particles which are also ill-defined in the Davies-Fulhng 
model in the presence of horizons. 
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Introduction 



When considering the fluxes emitted by a uniformly accelerated mirror described by the 
Davies-FuUing model, one encounters a paradoxical situation: when working in the initial 
vacuum, the local flux of energy vanishes whereas the Bogoliubov coefficients encoding pair 
creation do not. Moreover, (N^), the mean number of particles with frequency u emitted 
to J'^ diverges [|, |^, 0, and so does the corresponding total energy {H) = J^du Lj{Nt^). 

It should be pointed that similar properties are shared by all uniformly accelerated 
systems. They are indeed found (in a slightly different form) in the case of a uniformly 
accelerated classical charge @] , an accelerated two-level atom [§, , and for accelerated 
black holes |^ . The vanishing of the energy flux and the divergence of the total energy 
are both unavoidable when considering uniformly accelerated systems whose coupling to 
the radiation field is constant. The vanishing of the flux follows from the facts that the 
orbits with uniform acceleration are generated by the boost operator of the Lorentz group 
and that the Minkowski vacuum is a null eigen-state of this operator. Hence stationarity 
is built in and this implies a vanishing flux. (When applied to accelerated atoms this is 
known as Grove's theorem |Ty, [T^, 0-) The divergence of the total energy emitted 
follows from the singular behavior of the "Rindler" modes (i.e., the eigen-modes of the 
boost operator) on the horizons associated to a uniformly accelerated trajectory, see 
Appendix C of [|^ for a general proof. 

In order to obtain regular expressions for the flux and the energy emitted, one needs 
either to regularize the trajectory by decreasing the acceleration at asymptotically large 
times, or, what we shall do, switch off the coupling between the accelerated system (here 
the mirror) and the radiation field. This procedure was already applied to an accelerated 
two-level atom in and regular expressions have been found for the local flux and the 
total energy emitted. In the present paper, it is our intention to apply a similar treatment 
to an accelerated mirror. However, one immediately encounters a problem : this analysis 
cannot be performed within the framework of the original Davies-FuUing model because, 
by construction, the reflection on the mirror is total. Therefore, we first introduce an 
alternative model based on an action principle which, on one hand, reproduces the results 
of the Davies-FuUing model when the coupling between the mirror and the radiation is 
constant and, on the other hand, allows to switch on and off the interactions. 

When using this model to describe the scattering by an accelerated mirror, we obtain 
regular expressions for the transition amplitudes in the place of the singular Bogoliubov 
coefficients obtained with the Davies-FuUing model. Both local quantities, such as the 
mean flux (Tyy), and global ones, such as the mean energy (H) and density of particles 
(Ni^), are now well defined. Moreover, as long as the coupling is constant, we recover the 
fact that the energy flux vanishes. However, it is now preceded and followed by transient 
effects associated with the switching on and off. Because of the ever increasing Doppler 
effect associated with constant acceleration, these effects are exponentially blue-shifted 
(in terms of the proper time of the mirror). Therefore, in order to get a finite energy, the 
rate of switching off the interaction must be faster than the growth of the Doppler effect (a 
condition also found in [^]). If this condition is not fulfilled, the mirror emits an infinite 
energy, thereby recovering ill-defined results as those obtained with the Davies-FuUing 
model. 

To further clarify the physics into play in the scattering by a mirror of acceleration a, 
we compute the transition amplitudes governing pair creation in a "mixed" representation. 
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By mixed we mean that in each pair, one quantum is characterized by cu, a Minkowski 
frequency, whereas its partner is characterized by A, a Rindler frequency, the eigenvalue of 
the energy with respect to the proper time of the mirror. This representation is useful for 
the following reasons. First, when A = AM, the scattering amplitudes closely correspond 
to the transition amplitudes governing the absorption and emission of photons by an 
accelerated two-level atom whose energy gap is AM. It is quite nice to see how, for 
each Minkowski frequency uj, the "exciton" of the atom is replaced by the emission of the 
partner of the Minkowski quantum. Secondly, the range of Minkowski frequencies uj which 
participate to the processes when the interactions last a proper time lapse of T is finite 
and grows with ae"^. Hence the physics into play is a succession of pair creation acts 
with Doppler shifted Minkowski frequencies. However, this cannot be seen by considering 
only the expectation value of the flux, i.e. the one-point function of the stress tensor, 
because these successive individual effects perfectly interfere with each order since, in the 
absence of recoils, the orbit is generated by a Killing field. In order to reveal these local 
acts, one should either consider the correlations amongst emitted quanta by computing 
the two-point function of the energy flux or enlarge the dynamics so as to take into 
account the recoil effects. In a next article ||15[ we analyze these quantum correlations, 
and in a forthcoming paper we shall analyze recoil effects. 



1 The Davies-Fulling model and its extensions 

In this section, we first analyze the case of asymptotically inertial trajectories (Sec. 1.1). 
We introduce notations which also apply to the cases of trajectories which enter or leave 
space-time through null infinities (Sec. 1.2), cases which turn out to be more delicate to 
analyze. In Sec. 1.3 we introduce the self-interacting model. Sec. 1.4 is devoted to the 
study of energy fluxes. 

1.1 Asymptotically inertial trajectories 

In the Davies-Fulling model |]l|, one studies the scattering of a massless scalar field in- 
duced by imposing a Dirichlet boundary condition along a time-like trajectory in 1 + 1 
dimensions. The evolution of the field operator is governed by the d'Alembert equation 

{dl-dlmt,z) = ^, (1) 

together with the reflection condition along the trajectory of the mirror z = Zci{t) 

^{t,Zd{t))={] . (2) 

Since we work in 1 + 1 dimensions with a massless field, it is particularly useful to work 
in the light-like coordinates U^V = t ^ z. Then, Eq.(|l|) becomes dudy^iU, V) = 0, and 
its general solution is a sum of a function of U plus a function of V. Since the trajectory 
of the mirror is given once for all, the recoil effects induced by the scattering of quanta 
are neglected. 

In this sub-section, in addition to the condition that the mirror trajectory be time-like, 
we consider only asymptotically inertial trajectories, i.e. trajectories which originate from 
the time-like past infinity i' and which end in (see Fig. 1). Then, since the reflection 
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Figure 1: In this Penrose diagram, the sohd line is a time-hke trajectory going from to i+. 
The dashed hnes represent incoming V configurations which give rise to the production of a pair 
of outgoing U quanta (right movers). 

is total, the configurations emerging from J7^, tlie riglit part of J'~ , are completely 
decoupled from those emerging from J7^. Therefore, one can analyze what happens on 
each side separately. On the right hand side, all left movers are scattered into right movers 
and sent toward J7^. In second quantization, when the trajectory is not inertial, this leads 
to the production of pairs formed with two right movers. Similarly, on the left hand side, 
one studies the scattering from to J'j^. Since the expressions governing the scattering 
on the left are obtained from those on the right by exchanging U and V, we will restrict 
ourselves to the analysis of the scattering from to J7^. 

To analyze the scattering in second quantization, one first needs to identify the m-basis 
of modes which are defined before the scattering occurs. On J'^, the usual eigenmodes 
of Minkowski energy idt = uj > can be used since JT^ is a Cauchy surface for the left 
movers when the mirror emerges from i~ (this is no longer the case for trajectories which 
emerge from J7^, see Sec. 1.2). On J'j^ the Minkowski modes are given by 

^yM^v,u = = ^= , (3) 

where the upper index V, in means that the mode is left-moving and defined on the 
initial Cauchy surface . We have introduced the index V in order to be able to deal 
with partially reflecting mirrors for which left and right-movers should be simultaneously 
considered. The norm of v^^'*" is determined by the usual Klein-Gordon scalar product. 
On J7^, one has 

r+oo <_> 



where f dy g = fdyg - gdyf. 

For finite values of U, on the right hand side of the mirror, i.e. for V > Vci{U), the 
m-mode v^^'*", the solution of Eq.(|^) which has Eq.(^ as initial Cauchy data, is 

p-iujV „-iujV^i{U) 

^:^-iv,u) = ^-^^. (5) 
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To analyze its final frequency content, it should be Fourier decomposed on J^. In total 
analogy with what we have on J7^, on the positive frequency modes are 

^^r(f/,\/ = +oo) = -==. (6) 

Since they are complete, on JT"^, the in-modes can be written as 

VAnlTT t r I \ I J ' I UV * U.out nUV * U.out *\ 

[U, V = +00) = du [a^,^ - (3^,^ j . (7) 

When evaluated on J^^ the overlaps are given by 

r+00 



To interpret the scattering in terms of particle creation, one should decompose the 
Heisenberg field operator $ both in the in and owt-bases. When working with a complex 
field, annihilation m-operators of particles and anti-particles are defined by 

ar-^" = (v^r'^", <^>) , &r'^" = (v^r'",<f^) • (lo) 

Because the in-modes v^^'*" form an orthogonal and complete basis, these operators satisfy 
the canonical commutators when the field operator satisfies the equal-time commutation 
relation [^{t, z),dt^^{t, z')] = iS{z - z'). The in-vacuum state is defined, as usual, by 
a^'*"|0) = &^'*'^|0) = 0. Similarly, the out-operators are defined with the owt-modes V'J^'""*- 
Since we are dealing with a linear theory without sources, the overlaps a and (3 of 
Eqs.(^ and define the Bogoliubov coefficients relating the initial and final operators 
ai^,b^. Therefore, these overlaps determine the expectation values (as well as the non- 
diagonal matrix elements) of all operators built with $. For instance, when the initial 
state is vacuum, the mean number of right-moving particles of energy u' received on J'^ 
is 

poo 

(iv5) = .n(o|a^r*t^^r*lo),„= / du . (ii) 



— CO 

oo 



1.2 Non- asymptotically inertial trajectories 

When the trajectory does not end on (or does not begin from the strict decoupling 
between left and right movers is no longer realized even when the refiection on the mirror 
is total. Consider for instance the trajectory 

V,i{U) = -K-^e-^"". (12) 

The mirror goes from i~ to V = Q on . Depending on the sign of V ^ left-movers 
emerging from are either refiected into right-movers for \^ < or end as left-movers 
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Figure 2: In this Penrose diagram we represent the trajectory defined in Eq.(12). This trajectory 
has been often considered (see e.g. |l^) because of its analogy with the Hawking effect. In this 
case, is a Cauchy surface whereas is not. The portion of with V > Q plays the role of 
the future horizon of the black hole. The dashed lines are incoming left-movers. One sees that for 
V < the quanta are reflected, giving rise to right-movers for all values of U . On the contrary, 
for > 0, the incoming quanta do not reach the trajectory and end up in . The question 
mark is there to raise the reader's attention on the issue of the choice of the appropriate basis of 
out-modes to decompose the field configurations when the mirror crosses . 

on for V > (see Fig. 2). Thus on , the image of v^J^'*" of Eq.(^ now contains 
two pieces 

^v,n^u, V) = e{v)^ - . (13) 

This mode is singular at ^ = where the mirror hits J'j^. On the other side of the mirror, 
the U in-modes emerging from jT^f are fully reflected into left-movers but they are also 
singular on J'^. Hence both sets of m-modes are singular on JT""*", at = 0. 

This raises an interesting question: Given that the mirror trajectory ends on J'^ , 
which is part of the final Cauchy surface, what is the appropriate set of owt-modes to 
describe the scattered field configurations ? 

The procedure we shall follow is to decouple asymptotically the radiation field from 
the mirror, i.e. to make the mirror asymptotically transparent. In this case, the free 
Minkowski modes e"*"^'^/ y/A-Kuj' and e-'^^'^/VA TCLu' still form a complete owt-basis. We 
shall not adopt the other choice which consists in working with out-modes defined on 
either side of the mirror on J'j^ . These modes are so singular that their (Minkowski) 
energy content is not defined. Nevertheless, when working in a state specified on J'~ , 
the expectation values of local operators whose support is 7^ is independent of the 
out-basis one chooses. The out-basis is necessary only for computing global quantities 
such as the total energy {H) = J^du lj{N^). 

When adopting the asymptotic decoupling hypothesis, the image of v^^'*" on J'^ U 
{iJiiy > 0)} can be decomposed as 



V,in — f rli f n,^^ * ,n^'°^* /qJV* j,out*\ 

^ - du \^a^,^ ip^, - p^,^ ip^, j , (i4j 
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where v'^f"* are the usual Minkowski modes, as in Eq.(^. The Bogohubov coefficients 
a, (3 are now 2x2 matrices. The discrete index j stands for f/, V and is summed over 
when repeated. These coefficients are still given by the Klein-Gordon scalar product as 
in Eqs.(|D and (P), since the out-basis is composed of usual Minkowski modes. When the 
trajectory is asymptotically inertial, we recover what happens in the right hand side of 
the mirror for a"^ = a^^ , = and on the left hand side for a^^ = a^^ , = (3^^ . 
In addition, one has a^^ = (3^^ = a^^ = (3^^ = 0. Finally, we mention that a similar 



decomposition to Eq.(14) holds on each side of the mirror when the mirror travels from 



J' to JT""*", as it is the case for a uniformly accelerated system, see Sec. 2. 
1.3 Partially transmitting mirrors 

In preparation for subsequent analysis, we now present how to study partially transmitting 
mirrors. In this case, one should also consider simultaneously U and V modes. Indeed, 
when the trajectory is asymptotically inertial, an incoming U mode is partially scattered 
into an outgoing V mode and partially transmitted as an outgoing U mode. Hence, when 
the mirror is not inertial, a proper description of the Bogoliubov coefficients requires to 
consider 2x2 matrices a, (3 which mix U and V modes. 

There are two different ways to describe partially transmitting mirrors. First, one can 
choose from the outset the transmission coefficient (expressed in the rest frame of the 
mirror) and deduce from it the Bogoliubov coefficients, see Sec. II. B in |^. We shall not 



follow this method since it does not allow to switch off the coupling to the radiation field. 

The other method is based on self-interactions described by an action. The principle 
usefulness of this model is to allow to switch on and off the coupling of the radiation with 
the mirror. We will see in the next sections that this is necessary to obtain well defined 
transition energy fluxes for a uniformly accelerated mirror. In the following, we shall use 
only this model, see Sec. Ill in |]16| for more details. 

In this model, the scattering on the mirror is governed by an action whose density is 
localized on the mirror trajectory a;^;(r), where r is the proper time. 



-'int 



J dTH,^t{r) = -90 J drg{r) J (fx ^^(x^ - <,(r)) z) , z)] . (15) 



go is the coupling constant. The real function g{T) controls the time dependence of the 
interaction: When the coupling lasts a proper time lapse equal to 2T, g{T) is normalized 
by J drgij) = 2T. To preserve the linearity of the scattering, J-' must be a quadratic 
form of the field $ and to have a well-defined Hamiltonian, it should be hermitian. The 

various possibilities with the lowest number of derivatives are J-'o = JFi = ^Ud^^ 
and J^2 = dr^^dr^. 

In the interacting picture, the charged field evolves freely. It can thus be decomposed 

as 

V) = {ale-^ + ale-^ + h^^e^^^ + h^e^-^) . (16) 

Jo ^/^TlUJ 

The annihilation and creation operators of left and right-moving particles (and anti- 
particles) are constant and obey the usual commutation relations 

k, = 5''5ico - CO') , [bl, bil] = 5^^{co - CO') . (17) 



7 



All other commutators vanish. In the interacting picture, the states evolve through the 
action of a time-ordered operator Te'^'"*. When the initial state is vacuum, up to second 
order in qq, the state on J'~^ is 

Te^^-|0) = \0)+tL,nt\0)+^^^^f^\0) + \D) . (18) 
The ket \D) contains terms arising from time-ordering. None of these terms contribute to 



the total energy emitted (see [16| for a detailed analysis). Hence we drop \D) from now 
on. 

The relationship between this model and the original Davies-Fulling model can be 
made explicit by considering the case where J-' = ^Hdr^ and (?(r) = 1 (see JlBl). In this 



case, whatever the mirror trajectory is, the first order transition amplitudes are related 
to overlaps a^^, , /Jj'^, entering in Eq. ([I^) in the following way 

^L^ * = (0|ar e^^-* ^r^^lO), = Siuj - u') - zg^ (19a) 

Srj* ^ (0|e-^- aZ^bi;\0) = -.^o/^rj (19b) 

* ^ (0|ar e^^- a^J\0)^ = -ig^ al^^, * (19c) 

B::3*^{0\e-''-'aZ^b''J\0) = ^goP^3* , (19d) 



where the subscript ()c means that only the connected graphs are kept. In Eqs.([T9|b) and 
(p!9|d), one sees clearly the link between the (3 coefficients and pair creation amplitudes. 

When using these amplitudes to compute energy fluxes one encounters severe IR di- 
vergences due to the massless character of $. These divergences can be eliminated by 
considering the Hamiltonian with one more derivative 

^[$t^ $] = ^^(a^$ta^$ + a^$a,$t) . (20) 

ar ar 

The two terms within the parentheses mean that the interaction is symmetrical under 
charge conjugation. This implies that the transition amplitudes will be invariant under 
the exchange of a and b operators. 

It should be stressed that the self-interacting model can handle without modification 
in the cases when the trajectory enters and/or leaves space through null infinities. In these 
cases, when computing perturbatively transition amplitudes, one automatically adopts the 
convention of using asymptotically free modes. Indeed, the interacting picture is based on 
the assumption that the interaction is switched on and off asymptotically. This remarks 
reinforces the well-founded character of the choice adopted in the former sub-section to 
use, for the owt-basis, free Minkowski modes on JT""*" . 



1.4 Energy fluxes 

In this subsection, we compute physical observables such as the number of emitted parti- 
cles, the energy and its fiuxes, to second order in go and when the initial state is vacuum. 
We study only the left-moving quanta emitted toward J7^: all the results for the right- 
movers can be obtained by exchanging U and V. 

The mean number of V particles of energy u is particularly simple because only the 
second term of Eq.(|T8|) contributes. One obtains 

poo 

{N3) ^ (0|e-^- ar tar e^"-'|0) = du' {\B^3\' + \b33\') , (21) 
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in the place of Eq.(]Tl|) since the partner of can be either a U or a.V quantum. Then 
the (subtracted) integrated energy is, as usual, 



oo 



Hi) = 2 duj u{N^) , (22) 



where the factor of 2 stands for particles + antiparticles. 

One can also compute the local flux of energy. The corresponding Hermitian operator 
is Tvv = dv^^dy^ + dv^dv^^ ■ Its expectation value is given by 

{TvviV)) = (0|e-^^-* Tvv e*^-'|0) - {0\Tvv\0) 

= {nv) + {T¥v). (23) 

where 

{Tvv) = {^\LintTvvLint\0)c 

- 2J2 [ r i^dJ^e-'^----^- ( Hdk Bll -E^ , (24) 



j=uy 



and 



(T^V) = -2 [ \m{{Q\TvvLintm + Re {(0|TvyL,„tL,„t|0) J] 



We have subtracted the average value of Tw in the vacuum in order to remove the zero 
point energy. In Eq.(p3|), we have introduced {Tvv) which determines the integrated 
(positive) energy {H^) of Eq.(^) and {Tvv) which integrates to 0. Note that the linear 
term in qq in the first line of Eq . (^5]) reappears in the second through the definition of 
the A terms given by Eqs.(^9|). 

Besides these expressions based on the amplitudes A and B, one can also express 
{Tvv) by using the "scattered" Wightman function^. This function is defined by 

W{U, V- U', V) = (0|e-^^-* $t(f/, V) $(f/', V) e^^-*|0)^ . (26) 

To obtain the subtracted flux, one needs also the unperturbed Wightman function eval- 
uated in the vacuum 



W,ac{U,V;U',V') = (0|$t(f/,y) $([/', \/')|0) 

= -^(\n{V' -V -le) +\n{U' -U -ie)) . (27) 
47r 

■^This illustrates the fact that one does not need to choose an owt-basis when computing expectation 

values of local operators with initial states prepared on J'~ . Notice however that the subtraction in 

Eq.(p8|) implicitly reintroduces the notion of Minkowski vacuum on JT^"*" since the only singularity of 

Wyac{U, V; U', V) in Eq.(p7|) is the usual short distance one, independently of the presence of the mirror 

on ^7+. Thus, since subtracting the vacuum contribution on JT"^ is equivalent to subtracting that of 

the out- modes, the use of Eq.(p7|) implies that the mirror is no longer coupled to the radiation field at 

asymptotically late times. 
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In terms of these functions, the mean flux on Jj^ reads |T6| 



2 hm^ dvdv [WiU, V- [/', V) - W,ac{U, V- U', V')] . (2J 



Notice finally that this expression also applies to the Davies-FuUing model and leads to 
the well-known result 



(TwiV)) 



DF 



— lim hmdvdv'[HUciiV')-UciiV) 
2tt v'^v €~^o 



\n{V' - V 




dU, 



d 



dV 



-1/2- 



(29) 



From this equation one sees that the energy flux is local in that it contains at most three 
derivatives of UdiV) evaluated at the advanced time V. When considering the interacting 
model with g^r) ^ const., this local property will be lost. 



2 Uniformly accelerated mirrors 



Uniform acceleration means that 

^^ = -a^ = const. (30) 

In terms of Minkowski space-time coordinates, the trajectory reads t^—z^ = UV = —l/a?. 
In the sequel, we will consider a uniformly accelerated mirror living in the right Rindler 
wedge i?, i.e. its trajectory is UdiV) = —l/a^V with V running from to +oo, see Fig. 
3. 




Figure 3: In this Penrose diagram, we show a uniformly accelerated trajectory going from 
V — Q on Jj^ to J7 = on J^. The dashed lines represent the scattering of a pair of quanta 
(represented by localized wave-packets) emerging from ■ One particle is reflected into an 
outgoing V quantum for f7 < whereas the other ends as a left-mover on for [/ > 0. The 
dotted line represents a 1^-quantum which is not reflected by the mirror. 

The scattering associated with this trajectory leads to several difficulties when using 
the Davies- Fulling model. In Sec. 2.1, we list these difficulties. Then we will see how 
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they can be resolved by using our self-interacting model with a smooth switching on and 
off coupling. 



2.1 The difficulties using the Davies-FulUng model 

When considering the scattering by a uniformly accelerated mirror with the Davies- Fulling 
model, on the left side of the trajectory, the image on of the scattered in-modes are 

^l^^-iu^v) = e(f/)-= + e(\/)(-^-=) (31) 



^:''-{U,V) = e{-V)^. (32) 

As expected from Sec. 1.2, they are singular where the mirror enters (V^ = 0) and leaves 
{U = 0) the space-time. Their overlaps with plane waves (owt-modes) are 

= i^^r^ (33b) 

V /An y/ujuj' uj + uj' 

<ZZ*^ Ur\ = ^K,H^-^) (33c) 

Pi^'i^ =['^u' ^V^J j = — (33d) 

V / ira a 

where Ki{z) is a modified Bessel function (see Appendix A). 

When computing {N^} (given by J^du' HP^^H' + \P^FJ')), see Eqs.(0) and (|l])), the 
integral over u' diverges in the JR. Moreover, diverges when u = uj' . Similarly, a^,^ 
is ill-defined since the integral representation of the Bessel function in Eq.(^3]c) requires 
to contain a finite an positive real part, see Eq.(p3D. 

In addition to these problems in momentum space, when computing the space-time 
properties of the flux, one encounters the following properties. When plugging UdiV) = 
— 1/a^V for y > in Eq. ([29|) , one finds that {Tyy{V))^^ vanishes^. This is not in agree- 
ment with the non- vanishing character of the (3 since, on one hand, (H'^j) = j^duj uj{N^) 

and on the other hand, {H'^j) = f^^dV {Tyy). It is as if the created particles were car- 
rying no energy flux |l|, |T^, [l8l . 



-iuV 



To sum up, the difficulties of the Davies-FuUing model for uniformly accelerated tra- 
jectories are: 



unregulated and IR diverging overlaps, Eqs.( 
a diverging expression for (N^^), the mean number of particles created, 
a vanishing local flux although pair-creation transition amplitudes do not vanish. 



'^Generally, the trajectories which provide a vanishing V flux are of the form Uci{V) — ^y^^ 
Time-hkeness imposes AD > BC. If C = 0, we recover inertial trajectories. If C 7^ 0, we recover 



uniformly accelerated trajectories with a = C/y/ AD — BC. 
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when one considers the scattering by two uniformly accelerated mirrors with sym- 
metrical trajectories, i.e. which both obey UV = — the (unregulated) overlaps 
f3 also vanish, together with the local flux. From [|19|, one could infer that these 



settings form a perfect interferometer. This cannot be exactly the case since the two 
mirrors live in two causally uncorrelated regions. This issue will be fully discussed 
in a forthcoming paper ||15||. 



2.2 The switching function g 

To avoid the difficulties listed above, we shall use our self-interacting model, based on 



Eqs.(p!5|) and (^Ol) . We require that the switching function g{T) be continuous, differen- 
tiable and that it decrease sufficiently rapidly for large r. A choice we find very convenient 
and shall adopt is 



-2?7 cosh(ar) 



(34) 



where < r/ ^ 1 is a dimensionless parameter which controls the switching time. 

g{T) can be interpreted in two different ways: either as governing how the interactions 
with the mirror are turned on and off, or as a mathematical regulator which properly 
defines the transition amplitudes A'^j^ , Bfj^ . This second interpretation implicitly relies 
on the limit r] 0, see Appendix B. In the body of the text, we shall use the first 
(physical) interpretation of 5'(t). 




Figure 4: Here are presented two ^(t), both for a = 1. The sohd hne is with Inrj — —20 and 
the dashed line with hnj — —18. One sees that the size of the plateau is linear in In 77 whereas 
the slope is independent of rj. This remark will be crucial when considering the Rindler energy 
emitted by the mirror, see Sec. 3. 



The features of g{T) are the following (see Fig. 4) : 
a) a plateau of height 1 centered around r = and of width equal to 

2|ln(2r/)| 



2T = 27Tgx-. 



(35) 



b) slopes which are maximal and equal to a/e + 0{7f) for ar ~ ±(aT + In 2), 
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c) an exponentially decreasing tail. We shall see that this extremely rapid decreasing 
behavior sufficient for having a finite Minkowski fiux on . 

We now compute the Rindler and Minkowski Fourier components of g since they will 
reappear in the expressions of the transition amplitudes. They are given in terms of the 
modified Bessel functions: 

1 /"^"^ 1 
9x = — dT e-2^™^'^(-)e^^^ = —K,x/a{2v) (36) 
2n .1^ an 



9u 



1 /"OO 1 1 

= 7^ ^-viaV+l/aV)^^u.v ^ Ki{27] ^ 1 - tuj / ar]) . (37) 



In the UV, for A ^ a and u ^ a/rj, the Rindler and Minkowski components decrease 
respectively as 

\gx\ -(a/A)i/2g-.A/2a (38) 

a 



\9c^ 



UJ—>00 



-{a/oof/y/^e-V^ . (39) 



From Eq.(|38|), one sees that the UV behavior of the Rindler component is independent 
of 7]. This is a direct consequence of point b) which states that the maximal slope of the 
switching function is independent of rj when expressed in proper time r. From Eq.(|39|) 
one finds instead that the UV behavior of the Minkowski components is damped by the 
regulator rj. 



2.3 Regularized amplitudes and particle content 



Given g(T) of Eq. ([5^) , the transition amplitudes can be explicitly calculated. They are 
given in the Appendix B and, as expected, they are well-defined. Nevertheless, we will 
not work with these amplitudes characterized by Minkowski frequencies since they are 
not convenient to compute the expectation value of observables. Similarly, we will not 
work with the transition amplitudes with Rindler frequencies even though they are simply 
expressed in terms of gx of Eq. ([36| ) . 

It turns out that it is more convenient to express the fiuxes and the energy in terms 
of transition amplitudes containing one Minkowski and one Rindler frequency. More 
precisely, these amplitudes mix Minkowski and "Unruh" quanta. The Unruh modes (^■^ 
are linear combinations of positive frequency Minkowski modes and eigenmodes of Rindler 
energy A (see 0, |13[ and Appendix C). 



These "mixed" transition amplitudes^ are given by the matrix elements of the scat- 
tering operator with the Minkowski operator a^^ and the Unruh one d-^: 

Bl = {0\alP,e^^^-^\0) . (40) 



^These scattering amplitudes make contact with the transition amphtudes of a uniformly accelerated 
detector coupled to the radiation field <& ||, |l3[ ^ . For instance, the spontaneous emission amplitude 
of a two-level atom with an energy gap AM is equal to B^^\x=am , see Eq.(2.48) in 
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To first order in qq, using Eqs. (]I5|) , (|20|) , (p3|) and (|78D, we get 



= ^ y , _ f,../. (1 - .c./ar,)-(^+-Va)/2 ^,^^,,^(2^/13-7^) . (42) 

Moreover, as shown in the Appendix C, U and V Unruh modes coincide when evaluated 
along the trajectory. Therefore B^^ = B^^ and similarly for A. 

Using these amplitudes, the mean Minkowski energy emitted to J'^ can be written 

{Hl,)=A duu d\\BZ^\^= du dXhMiuj,X;v) ■ (43) 

^0 J -oo Jo J -oo 

The number 4 before the first integral means that processes involving U ot V particles 
and U OT V anti-particles equally contribute to the mean Minkowski energy emitted to 
(or J^). 

The exact computation of {H'^) will be performed in the next section. In the following, 
we shall study hM{^^, A, 77) = [-B^a'I^ to give a qualitative understanding of {H\j). The 
usefulness of the mixed representation is that the behavior of hM^uj, A; 77) in the (a;. A, rj) 
space is quite easy to explain. The starting point is that the norm of -B^^ decreases as 
g-27r|A|/a ^ gxpected from the correspondence mentioned in the last footnote. 

Hence, the relevant range of A is centered around and of extension a few a's. When A 
belongs to this interval, the following analysis applies. 

First, 7] acts as a regulator: the Minkowski frequencies which contribute to Hm belong 
to the interval 

aV<^< iah , (44) 

where ^ is a numerical factor. Its value is ~ 0.50 when one uses the mid height criterion: 
hM^uj = C,a/ri, A; ri)/hM{yJ = a, A; 77) = 1/2 with A belonging to the relevant interval. 

Secondly, within the range given in Eq.(^), Hm hardly depends on u, as shown in 
Fig.5.|^ Hence, for any given value of rj, one can first trivially perform the integral over 
Lj from ar] to ^a/rj. The value of this integral is given by hM^uj = a. A; 77) x ^a/rj, since 
?7 < 1. 

Thirdly, the height of the plateau hardly depends on rj, see Fig. 6. This can be under- 
stood from Eq.(p6|): when Eq.(|4^) is satisfied and when 77 <C 1, B^^ is independent of rj. 
Hence, one can take the limit 77 — to estimate how depends on A. In this limit, the 
scattering amplitudes obey 

B^x -m ^^/,_ ,,,, 7a. , (45) 



^In this, we recover what was found for an accelerated detector, see Sec. 2.4 of [jlsj. In that case, 
the Minkowski frequencies which contribute to the processes are the Doppler frequencies which resonate 
with the energy gap AM when the interactions arc turned on, i.e. those which satisfy lo = e^'^AAf for 
—T < T < T. Moreover, within that range, the transition amphtudes do not depend on to and the hmit 



?/ ^ can be used to estimate the amphtude, as in Eq.(45) 
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V -12 

-10 



Figure 5: Here we show hM{i^,^ = 0.15; 77) in terms of Ino; and In 77. One sees clearly that 
the surface exhibits a "plateau" of constant height which is limited by the lines In w = ± In 77 . 
/ijv/('^, A; 77) is given in arbitrary units and a — 1. 




Figure 6: Here we show hM{i^ — 1, A; 77) in terms of A and In 77. One clearly sees that the energy 
density Hm is independent of rj. This comes from that the fact that we compute it "within the 
plateau" , i.e. for ai] <ti lo = a ^ a/7;. Again, hMi^^, A; 77) is given in arbitrary values and a = 1. 



where 7]^ = {ip^, 0\) is given in Eq. (|77|) of Appendix C. In this case, we get 



hM^uj = a, X;r]) 



Auj\B. 
9l 



X 



27ia \ e^^l'' - e-'^Va 



(46) 



thereby recovering the above mentioned Boltzmann expansion law for |A| ^ a. Thus the 
mean energy can be approximated by 



e- 

V 



dX hMijjj = a, X; rj) 



67777 



(47) 



Although the numerical factor is not exactly determined (because of the ambiguity of 
defining C,), will see in the next section that the factor QqQ^ and the scaling in l/rj 
correctly define the behavior of {H'^,j) when r/ ^ 1 (The exact value of is 3/8 instead of 
1/2.). 
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So far, by the introduction of the switching off function (Eq.(plD), we have solved the 
first two problems listed in Sec. 2.1. Indeed, the transition amplitudes (expressed in the 
Minkowski, Rindler or mixed representation) are now well-defined both in the IR and the 
UV, and the energy emitted is no more infinite. The third and last point, i.e. the issue 
concerning the local fluxes, is the subject of the next section. 

3 Local Minkowski and Rindler fluxes 

In the Davies- Fulling model, the flux is given by Eq.(^). It is a local function of the 
trajectory UdiV) and its derivatives expressed at the advanced time V . This result relies 
on the time independence of the coupling. Indeed, this feature no longer occurs when the 
coupling to the radiation field is switched on and off. 

To compute (Tyy) of Eq.(p3D we first put together the two terms quadratic in qq. 
This is appropriate when computing local properties in space-time because it leads to a 
simplification since this gives rise to a commutator which is local. Then the flux reads 

{Tyviy)) = -2Im{(0|TyvIv.nt|0)} + Re((0|L,„4Tvy,L,„t]|0)) (48) 

= {Tyy{V)), + {Tyy{V)),. (49) 

(Note that all disconnected diagrams automatically cancel in this expression.) Both terms 
are governed by the second derivative of the Wightman function, Eq. (pTj) , 

3,3„W,,JV-V') = -l^^-^—-^^. (50) 
Using this function, they read 

{TvviV)), = -f^Im / dr g{T)V^,{T)-—-——--\ (51) 



Uci{r)Uci{r') 



+ 




where a dot designates a derivation with respect to the proper time. The function dy5 
comes from the commutator: 

[dy^\ dy.^] = [dy<l>, 9v"$^] = '-dy6{V - V) . (53) 

To evaluate the integrals, it is appropriate to use the dummy variable V = Vci{t) and 
to define a new function 

GiV) = VM[V])9ir[V]), (54) 

which can be interpreted as the effective coupling constant when using V as the time. 
Using this function, one can evaluate Eqs. (|51|) and ([5^ ) by integrating by part until the 
exponent of the pole is unity. All boundary contributions vanish if g{T) decreases faster 
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than e""'"^', a condition satisfied by the switching function we chose in Eg. (|3^ . If g(r) 
decreases slower than e"'^'""', the expectation value of Tyy is ill-defined. Hence it appears 
that the condition (7(r)e"''^' — > for r ±00 is a necessary condition for having well- 
defined Minkowski expressions. 

Concerning Eg. (|5lD , after three integration by parts, the last integration is trivially 
performed by using Im{(a; — = 7r6{x). Concerning Eg. (p2D , the two terms within 

the parentheses are egual. In order to compute this expression, one first perform the 
integral over r by using the function dyS. Then, as for Eg . (|5ll) , one integrates by parts 
until one gets a first order pole. 

Grouping the results for Egs.(0) and (|52D, one obtains, for V > 0, 

{TvviV>0)) = -^{godl.G-gliGdt.G + 2dvGdlG)) (55) 
- {idl.Gr + dy[...]) . (56) 



12n 

For V < 0, one gets {Tyv) = as expected since the V < part of J'j^ is causally 
disconnected to the mirror trajectory^. 

In Eg.(^), we have separated the flux into two parts, a sguare term which will lead 
to a positive Minkowski energy (H'^.j) = J^^dV {Tyy) and a total derivative which does 
not contribute to it (note the similarity with {Tyy) and {Tyy) of Eg. (|23|) ). When using 
the coupling function of Eg.(|3^ and Egs. ([56|) and (pSf), one can obtain an analytical 
expression for (-ff]^) 

(^m) = |^(8r7=^i^o(4r/) + (V + 2r/^)iri(4r^) 

-%7]^K2{^7]) - 3v^K:,{Ari) + ri^MAri)) (57) 

When taking the limit 77 — > one obtains 

^Hvyz'^l^l^om. (68) 

Hence, up to a numerical factor, we recover the result of Eg. (^7\) . It is interesting to see 
how the pathological features of constant coupling re-emerge when taking 77 — > 0. In 
this limit, the effective coupling constant of Eg. (^4]) obeys GiV) = aV. Hence Eg. (|55|) 
gives a vanishing flux whereas {H'^) clearly diverges (see Eg. (p8D ) . 

To complete the analysis of the transients, we now compute the Rindler flux ((T^^) = 
{dV/dv)^{TYv))) in terms of the Rindler advanced time v = ^\n{aV). This analysis 

''We would like to briefly comment on causality. When computing the flux in a causally disconnected 
point with respect to the trajectory, one must obviously find zero [^0|. This is trivially the case when 
one expresses, in Eqs.(p4[) and (|2^), the transition amplitudes as integrals over the proper time r and 
first performs the integral over lu. However, in the absence of a regulator, one looses causality when 
inverting the order of the integrations. This can be seen from the unregulated amplitudes where the fact 
that mirror was in the R or L quadrant is lost, see Eqs.([7l|) in Appendix B. The advantage of g{T) is to 



give rise to transition amplitudes, see Eqs.(|68|), wherein the prescription of the pole governed by 77 keeps 
control on causality. The same remark applies to the amplitudes in the mixed representation defined 
in Eqs.(^l|) and (^2|). When using them to compute the flux, causality is kept. Because of the analogy 
with the transition amplitudes of an accelerated detector, causality is also preserved by regularizing these 
amplitudes, see Eqs.(25),(26) and (63) of [§. 
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clearly establish that the Rindler energy carried by the transients effects is insensitive to 
the duration of the interaction and only depends on the rate of switching on and off the 
interactions. From Eq.(p5[), we get 



{Tvv{v)) 



12n 



127r 



- godvg + Qo {gdlg + 2{d^gf) 
godlg + gl {gd^g + ^d^gd^g) 



\2tx 



{a\d,gf^[dlgf)^d^\..\ . 



(59) 
(60) 



As for Minkowski energy (i?]^), Eq.(^) shows that one obtains also a positive Rindler 




0.5- 
0.4- 




CTo 



Figure 7: Here is plotted (T-ui,(u)) as a function of av = In(ay), for a = 1 in arbitrary units 
and for two different values of rj. The switching function has been taken for In?/ = —6 (plain 
curve) and Inr; = —8 (dashed curve). One can notice the two following properties. The flux is 
significantly non-zero within the transients only. The amplitude of these transients is independent 
of duration of the interactions governed by "q. 

energy {H^) = f^^dv (T„„). Using g given by Eg. (p^ , the Rindler energy is 

(K) = '-^v'K.iAv) (61) 

In Fig. 7, we have plotted the Rindler flux when the switching function g is given 
by Eq.(0). We previously noticed that the slope of g^r) does not depend on rj. As the 
slope determines the Fourier content of gx, one understands why we obtain a non- vanishing 
Rindler energy even in the limit rj —>■ 0. We could have chosen different switching functions 
g{T) such that their slope would tend to zero. In the limit, we would have found (Hn) = 0, 
as for a constant coupling. However, in this case, we would have necessarily obtained a 
diverging Minkowski energy (H^), since the condition (yf(r)e'^^'^' — >■ for r ±oo would 
not have been fulfilled. Hence {Hji) cannot be sent to if one requires to have a finite 
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Conclusions 



By considering the self-interacting model defined by Eg . ([15|) , with T given in Eq.(poD and 
(?(r) given in Eg. (jS^) , we have solved all the difficulties listed in Sec. 2.1: The transmission 
amplitudes are well defined and given in Egs. (^ID , (|42D and (|68|) , the mean energy flux is 
given in Eg.(|57D and the local flux in Eg.(^). All these guantities are regularized by the 
parameter r/ which controls the switching on and off of the coupling through Eg.(0). 

The important lesson which emerges from this analysis is the following: when express- 
ing (Tyy) in terms of A and B as in Egs.(El) and (ESf), the regulator r/ should be sent to 



after having integrated over /c, a; and uj' . In this, we recover what was found in when 
evaluating the energy density in the Rindler vacuum. If instead one first sends 77 — > 0, 
the unregulated expressions of the scattering amplitudes are so poorly defined that one 
even looses causality and crossing symmetry, see Appendix B. Therefore, in the presence 
of horizons, or more generally when the mirror enters or leaves space-time through null 
infinities, it is mandatory to consider the scattering amplitudes as distributions and not 
only as functions of freguencies belonging to [0, +00 [. 

In addition, by expressing the scattering amplitudes in the mixed representation, see 
Egs.(^) and (^21), we have made contact with the physics of a uniformly accelerated 
detector. Indeed, its absorption/emission transition amplitudes are given by the same 
functions as the scattering amplitudes in the mixed representation. This strict correspon- 
dence establishes that the physics of uniformly accelerated systems is dominated by the 
kinematics, namely (near) stationarity with respect to proper (Rindler) time and (near) 
singular behavior due the exponentially growing blue-shift effects associated with uniform 
acceleration. 

Appendix A : Bessel functions 

In this appendix we recall some features of the modified Bessel functions Ky{z\ where 
(z/, 2) e C (see p.374). They can be expressed by the following integral representation 



/■oo 

K,{z) = / rfte-"™^*^Wcosh(z/t), (63) 
Jo 

where | arg(z)| < 7r/2. For k & Af and u one has 

1 f) \ ^ 

-—] {z^'K^iz)} = e-'-'^z-'-'^K.^kiz) , (64) 
zoz J 



(~)'{^-AV(.)} 



e-"z-'"''K,Mz) . (65) 



We also recall the asymptotic behavior of the K^s for small and large arguments 

iro(^) '~°-ln(^) and K,{z) '^{-Y , for 3?(z/) > (66) 



whereas 



2 'z' 



K{z) '^^^ V^^"'' for all z/. (67) 
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Appendix B : Regularized transition amplitudes 



The aim of this appendix is to give the exact expressions of the regularized transition 
amphtudes in terms of Minkowski frequencies. The main virtue of the regulator 7] is to 
define them without ambiguity. The direct evaluation of Eqs.(|l9D with g{T) defined by 
Eq.(pl) gives 



A 



vv * . 



■■S{uj-uj') 



UJU' 



IT 



■K2{X) 



where X = 1r]\/\ — — uo') / ar] 



B, 



vv . 



^Wq a/I 



UOLjJ'\ T] 



:K,{X') 



A 



vu * 



TT a X'2 
where X' = 1r]\/\ — i{uj + uj')/arj 

Wo V\ 



\UJUJ 



71 a 
where Y - 

vu _ Wo 



\UJUJ' 



TT 



a 



KoiY) 

-- 2a/ {u/a + ir]){—uj' /a — irj) 
Ko{Y') 



where Y' = 2^/{uj/a + irj) [uo' / a — irf) . 



(68a) 



(68b) 



(68c) 



(68d) 



Two important remarks should be made. Because these amplitudes have been regularized, 
they possess analytical properties which guarantee that they obey crossing symmetry, that 
is 



A'^, : = -B' 



(69) 



see 



23 1 for exploiting this symmetry in studying accelerated detectors. Secondly, had the 
mirror followed the accelerated trajectory in the left quadrant rather than in the right one, 
the corresponding transition amplitudes would have been obtained by simply replacing t] 
by -7]. 

We wish also to stress that g{T) can be considered as a mathematical regulator which 
properly defines the transition amphtudes. Consider for instance B^^,. Using Eq. (|3^) , it 
is given by 



UJUJ' 



}_ i {{u + i'q)V - {u' - ir])/a^V) 



(70) 



One clearly sees that the integral is now well defined both for V ^ and V ^ oo. 

Finally, it is interesting to take the limit — > to see how one recovers the singular 
amplitudes that one would have obtained with a constant couphng. Using Eqs.(^), in 
the limit ?7 — 0, we get 



A 



vv * 



a/|cJZ7[ 



r/ >^ , Wo 

■d{uj -uj) + — a 

2n [UJ — uj') 



'\2 



B 



vv 



'271 



(a; + uj'Y 



(71a) 
(71b) 
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^vu . _ _^ Vk^ (71c) 



TT 



B-J-.l^^M2^). (71d) 

TT a a 

Although the choice of our Lagrangian, based on Eq.(^), has removed the IR diver- 
gences, the A terms are clearly ill-defined, as in the original Davies-Fulling model. More 
importantly, crossing symmetry and causality are both lost if one uses these unregulated 
amplitudes. This clearly establishes that the Bogoliubov coefficients must be conceived 
as distributions, or at least analytical functions of uj and tu', and not merely as functions 
defined from [0, +oo[. Thus the limit — should be performed only at the end of the 
calculation, after having performed all integrations. This is because the limit 77 — in 
general does not commute with these integrations. 



Appendix C : The Unruh modes 

By definition, the "Unruh" modes and (p^ possess the following properties: 

• they are made of positive Minkowski frequency modes only, whatever the sign of A 
be, 

• they are eigenfunctions of iaVdy (or —iaUdu) with eigenvalue A . 

They are thus well-adapted to study uniformly accelerated systems since they are eigen- 
modes of idr = A where r is the proper time calculated along the accelerated trajectory. 

Since the Unruh modes form a complete and orthonormal basis, one can define in 
a canonical way the corresponding annihilation and creation operators of particles and 
anti-particles a\,a\\b\ and b\^: 

a\ = $ ) , 61 = $t) , (72) 

where the subscript * stands as before for U and V. Hence, the scalar field can be 
decomposed as 



/+00 
., -00 ^ 



(73) 



i=U,V 



Note that the integrals over A cover the entire real axis. 

The Unruh modes are analytically expressed by the following expressions 

^v^y^ . l.n. (74) 

V^47rA(l - e-2-A/a) 

D -iuiV 

du^ iL^^ , (75) 



with 



iL = (^:,^r) (76) 



a \/2nau 



(7^)^: . (77) 



A sinh(7rA/a) 
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Notice that the regulator e in Eqs. ( |7^ and ([77| ) plays a role similar to rj in the text : e is 



needed to properly define the energy density in the Rindler vacuum p2 |. 

When considering U modes, we get similar expressions with 7^^ = 7]^ . Finally, when 
evaluated along the accelerated trajectories, within the right (R) or the left (L) quadrant, 
U and V Unruh modes coincide and are given by 



V = (r) = e-/a ^.y _ ^.u^rn t 



-iXr 
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